Optical characterization of ground states of polyacene 
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We investigate the ground-state properties of polyacene in terms of an extended Peierls-Hubbard 
Hamiltonian with particular emphasis on its structural instability of two types: double bonds in 
a cis pattern and those in a trans pattern. Calculating the polarized optical conductivity spectra 
within and beyond the Hartree-Fock scheme, we reveal a striking contrast between the cis and trans 
configurations. The two Peierls-distorted states are almost degenerate in their energetics but quite 
distinct in their optics. 

PACS numbers: 78.30.Jw, 78.20.Bh, 63.20.kd, 78.20.Ci 



I. INTRODUCTION 

There is a lot of enthusiasm for semiconducting or- 
ganic materials, which are available at low cost, easy 
to process, lightweight, and flexible. Especially for the 
last few decades, particular attention has been paid to 
polyacene oligomers, consisting of an aromatic linear ar- 
ray, due to their multifunctional electronic properties. 
Tetracene and pentacene, for instance, serve as light- 
emitting devices 1 ^ and field-effect transistors^^ The 
recent discovery of graphene 6 and the success of pattern- 
ing it into a narrow ribbon 7 -^^ have stimulated a renewed 
interesti - in polyacene as a basic building unit of them. 

Prior to the advent of the modern microelectronics, 
not a few scientists were already interested in the energy 
structure of polyacene from a theoretical point of view. 
Early investigations, initiated by Salem and Longuet- 
Higginsjii focused on the ground state of polyacene — 
whether and how the Peierls distortion occurs. The dis- 
cussion was more and more activated with the fabrication 
of highly conductive irans-polyacetylene?! 2 . There may 
be a structural instability in polyaceneji 1 - as well as in 
polyacetylene,— on a qualitatively distinct mechanism, 14 
however. Although much effort has so far been devoted 
to predicting which structure is energetically preferable 
in polyacene — cis-distorted, irons-distorted, or uniform 
(cf. Fig. [T|), this long-standing problem is not yet fully 
settled. Every structure was actually nominated for the 
most stable in the literature^ The controversy suggests 
that all these structures should closely compete with each 
other. 

First of all the cis- and trans-distorted isomers 
of infinite polyacene have exactly the same energy 
within an electron-phonon model without any Coulomb 
interaction 1 ^, such as the Su-Schrieffer-Heeger (SSH) 
Hamiltonian for polyacetylene^ They remain degener- 
ate with each other unless we take account of Coulomb 
interactions beyond the Hartree-Fock (HF) scheme. Sec- 
ondly the Peierls instability in polyacene is conditional 
for both cis and trans forms ) 11 ' 18 where the electron- 
phonon coupling strength seriously affects the scenario. 
Thirdly electron correlations make the situation more 
difficult. 1 ? They enhance the Peierls instability on one 
hand, while they push the Mott-Hubbard antiferromag- 



netic insulator on the other hand, against the uniform 
aromatic configuration. 

In such circumstances, Ramasesha and coworker o 18 ' 19 
applied modern numerical tools to polyacene modeled 
on correlated Hamiltonians of the Hubbard and Pariser- 
Parr-Pople (PPP) types, which include the on-site repul- 
sion alone and the power-law-decaying long-ranged inter- 
action whole, respectively, but both take no account of 
spontaneous phonons. They claim that the cis distortion 
is lower in energy than the trans one, both of which may, 
however, be replaced by a spin density wave (SDW) on 
the uniform lattice with increasing Coulomb interactions. 
While other author s 20 ' 21 ! 22 predict a fully delocalized 
nonalternating nature of polyacene employing standard 
and developed density-functional methods, the energy 
difference among such competing states must be rather 
small 1 - in any case. Then distorted and undistorted re- 
gions may coexist in actual polyacene compounds due to 
inevitable local defects and/or possible thermal excita- 
tions. We are now led to take an interest in character- 
izing these likely states rather than ranking them. Even 
if there is a slight energy difference between the cis- and 
trans-distorted structures, how can we tell the lower from 
the higher in practice? We suggest the polarized optical 
spectroscopy for polyacene in answer to this question. 

It deserves special mention in the context of our in- 
terest that Sony and Shukla 2 ^ have recently calculated 
the optical absorption spectra of oligoacenes on a large 
scale within correlated Hamiltonians of the PPP type. 
Starting from the restricted HF state, which is an undis- 
torted paramagnetic metal, they investigate in detail 
configuration interactions (Cis) between its particle-hole 
excitations of various levels — single, double, and more. 
In order to describe the effects of electron correlations 
on the optical properties of aromatic molecules as pre- 
cisely as possible, their argument is restricted to small 
oligomers up to seven benzene rings and any possibil- 
ity of lattice distortion is ruled out. Other groups of 
theoretician a 16 ' 24 i 25 i 26 i 27 ' 28 i 29 investigate an infinite acene 
chain within uncorrelated models of the SSH type in an 
attempt to find nonlinear excitations such as solitons 
and polarons, where any scenario is constructed on a 
Peierls-distorted structure without any concern for its 
configuration — cis or trans — because of the degenerate 
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energetics. Thus and thus, there is, to our knowledge, 
no guiding calculation of the optical properties of poly- 
acene under both electron-electron and electron-lattice 
interactions adequately existent. 



II. MODELING AND FORMULATION 

We aim at revealing the generic behavior of poly- 
acene rather than listing individual features of small 
oligoacenes. In order to simulate an infinite chain effi- 
ciently, we adopt the periodic boundary condition and 
restrict Coulomb interactions to a certain range. Srini- 
vasan et alJ£ investigated the ground-state properties 
of polyacene applying a projector Monte Carlo method 
to the Hubbard Hamiltonian with periodic boundaries, 
while Raghu et aL— carried out similar calculations ap- 
plying a density-matrix renormalization-group method 
to the PPP Hamiltonian with open boundaries. They 
declare that the long-range nature of electron-electron 
interactions does not qualitatively affect their essential 
findings such as the conditional Peierls instability and 
the cis distortion favored over the trans one. Sony and 
Shukla 23 further demonstrate that even the Huckel mod- 
eling succeeds in reproducing some aspects of the linear 
optical spectra of fully correlated oligoacenes modeled on 
the PPP Hamiltonian. The low-lying optical absorptions 
are well describable with both Huckel and PPP models, 
though the Huckel scheme considerably underestimates 
the optical gap in general. 

Thus convinced, we depict polyacene in terms of an 
extended Peierls-Hubbard Hamiltonian: 
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where cj.j a and q :j)S (cj. ■ s c; : j :S = n; y -, s ) create and an- 
nihilate, respectively, a it electron of spin s =|,|= ± 



at site j on chain I, while ri-j denotes the bond dis- 
tortion caused by the jth and [j + l)th carbon atoms 
on the Zth chain, on the description of polyacene as a 
couple of the trans isomers of polyacetylene with inter- 
chain bonding at every other site. The Coulomb inter- 
actions, ranging over neighboring sites, are modeled in 
6117a^ 1 s (Ref. where K is a 
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dielectric parameter, while aii and a± are the average 
lengths in A of neighboring C— C leg and rung bonds, 
respectively. Keeping the screened Coulomb parame- 
ters in mind, which were initiated by Chandross and 
Mazumdar— and successfully applied to oligoacenes by 
Sony and Shukla?^ we adopt U = 8.0 eV with k = 2.0. 
Though the Coulomb interactions in our use are not in- 
finitely ranged but cut off, such a parametrization must 
be suggestive and convincing under long-range correla- 
tions of merely moderate effect^ The intrachain and 
interchain electron hoppings are described by tn and 
tj_, respectively. We take 2.4 eV for (Ref. l23f ) and 
set i ± equal to 0.864* | (Refs. HlH^Ellii) , consider- 
ing that on ^ 1.4 A < a± ~ 1.5 A (Refs. EUS). a 
characterizes the electron-lattice coupling with K be- 
ing the c-bond elastic constant. We use, unless other- 
wise noted, a = 4.1eV/A and K = 15.5 eV/A 2 (Refs. 

[HiililiHiSliiii). 

In order to calculate the polarized optical conductivity 
spectra of polyacene, we define current operators along 
the long and short axes as 
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Since the charge-transfer excitation energy is of eV order, 
the system effectively lies in the ground state at room 
temperature. Then the real part of the optical conduc- 
tivity reads 

a l± (cj) = - Y^mUiM^SiEt -E - foj), (2.4) 



where \Ei) denotes a wave vector of the zth- lying state of 
energy The state vectors are calculated within and 
beyond the HF scheme, being generally defined as 
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where |-Eo)hf = IIrn^i a m|0) ^ s tne ground-state HF 
wave function with |0) being the true electron vacuum = 
and a\ n creating an electron in the mth HF orbital of b 
energy e m . The orbital label m is a function of mo- 
mentum k, band label A, and spin s, which are all good 
quantum numbers here. Any transition of finite momen- 
tum transfer is optically forbidden, which serves to re- 
duce the number of configurations to take into calcu- 
lation. Every excited state of the HF type is a sin- 
gle Slater determinant f 3 - 3 - where f{k,^,,v,s]i) = Sk^s.i: 
f '(hi, k 2 , Mi,M2, Vi,v 2 ,Si, S2;i) = 0, and thus E t = 
HF (E \H\E ) HF _ £m(fe,/i,s) + £m(fc,i/,s)- Those of the CI 
type consist of resonating Slater determinants* 3 - 4 - where 
the coefficients are determined so as to diagonalize the 
original Hamiltonian (|2.ip . Within the single-excitations 
CI (SCI) scheme, f(ki, k 2 , (ii, (J<2, vx, v%, si, s 2 ; i) remains 
vanishing and therefore, no excited-state Slater determi- 
nant mixes with the ground-state one: \Eq) = |i?o) HF . 
When we proceed to the single-double-excitations CI 
(SDCI) scheme, there occurs a significant correction to 
the ground-state energy as well as to every excited- 
state one. The CI method enables us to systematically 
investigate many-body effects beyond the HF approxi- 
mation in fairly large systems that we can hardly di- 
agonalize directly. It was indeed successfully applied 
to oligoacenes^ 3 - and related phenyl-based conjugated 
polymer a 35 ' 36 i 37 ' 38 i 39 ' 40 with its expansion truncated at 
varying level. 
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FIG. 1: (Color online) HF (dotted lines) and SDCI (solid 
lines) calculations of a ground-state phase diagram on the cu- 
ll plain at = 16 under the screened-Coulomb parametriza- 
tion. x and o indicate the two sets of model parameters 
considered here. We find a metallic state of aromatic (AM) 
configuration, an antiferromagnetic (AF) Mott insulator, and 
cis- and trans-distorted Peierls insulators, which read as in- 
phase (IP) and out-of-phase (OP) bond-order-wave (BOW) 
states, respectively. All the phases obtained are schemati- 
cally illustrated, where various segments qualitatively repre- 
sent the variation of local bond orders, while arrows in circles 
depict the alternation of local spin densities. IP BOW and OP 
BOW are degenerate with each other within the HF scheme, 
whereas IP BOW is slightly favored over OP BOW in the 
SDCI scheme. 



III. GROUND-STATE PHASE COMPETITION 

Even though we take particular interest in the optical 
features of polyacene as a Peierls insulator, it is still im- 
portant for us to have a bird's-eye view of its competing 
ground states. We calculate their energies at the HF and 
SDCI levels and visualize the ground-state phase com- 
petition in Fig. [TJ Besides Peierls-distorted structures, 
which are here referred to as in-phase (IP) and out-of- 
phase (OP) bond-order-wave (BOW) states, we find two 
undistorted structures, a metallic state of aromatic (AM) 
configuration and a Mott-insulating state of antiferro- 
magnetic (AF) configuration. AM is fully symmetric and 
has no band gap at the HF level. There opens a gap in the 
SDCI scheme. AF is more stabilized and wider gapped 
with increasing on-site Coulomb repulsion. There is an- 
other undistorted structure possible, which is character- 
ized as a charge density wave (CDW)^i but it is not sta- 
bilized into the ground state under the scrccncd-Coulomb 
parametrization. CDW is realized when we adopt a 
standard-Coulomb parametrization, 23 which is, however, 
inferior to reproduce experimental findingsj 35 i 36 ' 37 i 38 ' 39 i 40 



The HF energies of IP BOW and OP BOW are the 
same. IP BOW gains more correlation energy than OP 
BOW in the SDCI scheme. However, their energy differ- 
ence is small and decreases with increasing conjugation 
length. IP BOW and OP BOW are highly degenerate 
with each other in sufficiently long acene chains. They 
are stabilized conditionally, that is, depending on the 
electron-lattice coupling strength, against AM and AF 
under weak and strong electron-electron correlations, re- 
spectively. Moderate Coulomb interactions enhance the 
Peierls instability. The realistic parameters in our use, 
which are indicated with x in Fig. [TJ sit in close vicin- 
ity to a phase boundary. The ground state is a Mott 
insulator, but it closely competes with Peierls insulators. 
Many-body electron correlations seem to contribute to- 
ward a closer competition between them. They are very 
much likely to coexist in polyacene4£ Then, how can we 
distinguish between the two distorted structures? 
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FIG. 2: (Color online) HF, SCI, and SDCI calculations of the long-axis-polarized optical conductivity spectrum for IP BOW, 
OP BOW, and AF at N = 16 under the screened-Coulomb and moderate-coupling parameters. Every spectral line is Lorentzian 
broadened by a width of 0.06t n . 



IV. OPTICAL CONDUCTIVITY SPECTRA 

We calculate the polarized optical conductivity spectra 
of energetically competing IP BOW, OP BOW, and AF 
with varying number of contributive Slater determinants 
and compare them in Fig. O First of all we are im- 
pressed with an optical contrast between IP BOW and 
OP BOW, which is noticeable even in the HF scheme. 
They are degenerate with each other in their energetics 
but distinguishable from each other in their optics. Their 
most intense peaks appear far apart and in between is 
that of AF. The HF description of the relative intensity 
and position of each peak is much poorer than the CI 
findings, but it is useful enough to illuminate the optical 
nature of individual phases in a qualitative manner. 

Based on the HF energy scheme, Fig. [3] analyzes the 
polarized optical conductivity spectra parallel and per- 
pendicular to the conjugation direction. There are four 



molecular orbitals in each unit cell and they are molded 
into two valence and two conduction bands fulfilling the 
electron-hole symmetry. The dispersion relations of IP 
BOW and OP BOW are exactly the same within the HF 
scheme and remain alike even with many-body Coulomb 
correlations fully considered^ Their dipole transition 
matrices are also the same provided the excitation light 
is polarized in the rung direction, as is shown in Fig. 
[2c), where most of the spectral weight comes from 2-to- 
3 interband transitions. It is due to the parity-definite 
molecular orbital o 14 i 43 that both l-to-3 and 2-to-4 in- 
terband excitations make no contribution to the rung- 
direction optical conductivity. Under the reflection about 
the plane bisecting every rung bond, the valence and con- 
duction bands labeled 1 and 3 are of symmetric charac- 
ter, whereas those labeled 2 and 4 are of antisymmetric 
character, in both IP BOW and OP BOW. Any tran- 
sition with parity unchanged is optically forbidden. On 
the other hand, photoirradiation in the conjugation di- 
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FIG. 3: (Color online) HF calculations of the dispersion relations of the -zr-electron valence (A = 1, 2) and conduction (A = 3, 4) 
bands (a) and the "momentum-resolved" polarized optical conductivity parallel (b) and perpendicular (c) to the long axis for 
IP BOW, OP BOW, and AF at N = 16 under the screened-Coulomb and moderate-coupling parameters, where the area of 
each circle corresponds to the spectral weight. 



rection reveals a striking contrast between IP BOW and 
OP BOW, as is shown in Fig. [3^b) , where all the spectral 
weight comes from l-to-3 and 2-to-4 interband transitions 
in IP BOW, while little contribution from them in OP 
BOW. The conjugation-direction optical conductivity on 
an OP-BOW background arises mostly from 2-to-3 and 
l-to-4 excitations. 

The a\\ (u) spectra of IP BOW and OP BOW are thus 
distinguishable. Their main peaks sandwich the AF low- 
energy absorption bands (Fig. With Coulomb corre- 
lations fully taken into account, the most intense absorp- 
tion band is sharpened and grown up in general. Those 
of OP BOW and AF are close to each other, staying 
away from that of IP BOW. The screened-Coulomb and 
moderate-coupling parameters in our use stabilize AF 
slightly more than IP BOW. Band gaps are barometers 
of stabilization and therefore AF is gapped wider than 
IP BOW [Fig. E](a)]. Although the BOW Peierls gap is 
smaller than the AF Mott-Hubbard gap, the main ab- 



sorption peak in the AF spectrum appears much below 
that in the IP-BOW one. This trick is due to the optically 
forbidden lowest-lying Peierls-gap excitation through a 
long-axis-polarized photon on an IP-BOW background. 

The optical contrast between IP BOW and OP BOW 
is more and more accentuated with increasing coupling 
strength. If we reduce the elastic constant K to half 
the present value, as is indicated with o in Fig. [T] the 
electron-lattice coupling is effectively strengthened and 
any distorted structure is sufficiently stabilized against 
the Mott insulator. It may be the case, for instance, with 
organic molecular compounds of alternating tetrathiaful- 
valene and chloranil 44 and/or halogen- bridged transition- 
metal linear-chain complexesi 45 i 46 Figure 0] shows the 
same calculations as in Fig. E] but K = 7.75 eV/A . 
In the case of electron-lattice interactions predominat- 
ing over electron-electron correlations, both valence and 
conduction bands are much less dispersive in general and 
the long-axis-polarized dipole transitions between them 
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FIG. 4: (Color online) The same as Fig. [3] but the screened-Coulomb and strong-coupling parameters. 



are really eloquent of their background lattice distortion. 
In the conjugation direction, l-to-4 and 2-to-3 dipole ma- 
trix elements vanish on an IP-BOW background, while 
those of l-to-3 and 2-to-4 are negligibly small on an OP- 
BOW background. For light polarized in the rung direc- 
tion, IP BOW and OP BOW similarly behave and their 
optical features are essentially given by 2-to-3 interband 
transitions. 

We present in Fig. [5] the resultant polarized optical 
conductivity spectra. For light polarized in the conju- 
gation direction, IP BOW and OP BOW exhibit a sin- 
gle and two well-separate absorption bands, respectively. 
With configuration interactions fully taken into account, 
the single IP-BOW absorption band consists of an intense 
peak and its satellite in the high-energy side, whereas 
a couple of the OP-BOW absorption bands are both 
single-peaked, in a practical sense, due to the predom- 
inant zone-center excitations. Many-body electron cor- 
relations further have a significant effect on the relative 
intensity of the absorption bands. The spectral weight of 
the higher-energy OP-BOW absorption band turns out 



much smaller than that of the lower-energy one. For light 
polarized in the rung direction, IP BOW and OP BOW 
are degenerate with each other. Their spectra slightly 
differ beyond the HF scheme, but the difference is hardly 
recognizable in practice. The common a±(uj) spectrum 
is peaked similarly to the IP-BOW <7|| (w) spectrum, but 
high-energy absorptions are strongly suppressed in the 
rung direction. 



Finally we stress that the present findings are not indi- 
vidual features of small oligoacenes but symbolize poly- 
acene in the thermodynamic limit. All the features but 
irregular outlines peculiar to small clusters in Fig. [5] 
remain unchanged with further increasing conjugation 
length, as is demonstrated in Fig. [6] Indeed the absolute 
location and the relative intensity of every absorption 
depend on the number of benzene rings^ but they are 
almost converging at N = 16. 
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FIG. 5: (Color online) HF, SCI, and SDCI calculations of the polarized optical conductivity spectra parallel (||) and perpendic- 
ular (_L) to the long axis for IP BOW and OP BOW at N = 16 under the screened-Coulomb and strong-coupling parameters. 
Every spectral line is Lorentzian broadened by a width of 0.06t||. 



V. SUMMARY 



We have optically characterized competing ground 
states of polyacene and revealed a striking contrast be- 
tween the two highly-degenerate Peierls-distorted states 
in particular. The cis- and £rcms-distorted structures, 
which we refer to as IP BOW and OP BOW, are 
almost degenerate in their energetics but quite dis- 
tinct in their optics. The highest-occupied- molecular- 
orbital (HOMO)-to- lowest-unoccupied-molecular-orbital 
(LUMO) transition through a long-axis-polarized pho- 
ton is allowed in OP BOW but forbidden in IP BOW. 
There appear two well-separate absorption bands against 
a background of OP BOW and in between does a single 
absorption band with an IP-BOW background sit. Such 
distinct features may not be demonstrated as they are 
in polyacene of predominantly strong electron-electron 
correlation, but with growing electron-lattice coupling 



they become detectable literally. We may consider het- 
eroacenes such as paracyanogen^i as well as substituting 
hydrogen atoms in polyacene with larger molecules, in 
an attempt to tune the elastic properties and to realize 
a coupling-dominant situation. 

The realistic set of parameters is located in close vicin- 
ity to a phase boundary, where an antiferromagnetic 
Mott-insulating state, which we refer to as AF, is slightly 
lower in energy and gapped wider than the two Peierls- 
distorted states. The most intense AF absorption peak 
nevertheless appears far below the IP-BOW absorption 
band for light polarized in the conjugation direction. For 
light polarized in the rung direction, on the other hand, 
the common optical gap of IP BOW and OP BOW is 
naively smaller than that of AF. The relative location 
of their main absorption peaks remains qualitatively un- 
changed with varying configuration interactions, though 
the excitonic effects on the optical spectra are significant 
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FIG. 6: (Color online) The same as Fig. [5] but N = 128, where the SDCI scheme is never feasible. 



in general. 

There is another example 34 of optically characteriz- 
ing distinct ground states in competition. Some class 
of platinum-halide ladder compounds exhibits Peierls- 
distorted CDW ground states of the IP and OP types. 
According to their interchain valence arrangements, their 
optical conductivity spectra arc differently peaked. In 
this case, however, IP CDW and OP CDW not only look 
different in their energetics but also possibly present some 
contrast for the Raman spectroscopy. In the present 
case, IP BOW and OP BOW are so degenerate with 
each other that much effort has been devoted to solv- 
ing the problem of which is energetically preferable. 15 
Then the optical contrast between them will indeed come 
in useful for identifying them. Regular substitution of 



carbon atoms in polyacene, for instance, with nitrogen 
atoms, may lead to further stabilization of the Peierls- 
distorted structures against the aromatic one. IP BOW 
and OP BOW remain closely competing with each other 
in polypyridinopyridine4£ We hope our calculations will 
stimulate a renewed interest in polycyclic aromatic hy- 
drocarbons. 
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